Thermal correction to resistivity in dilute 2D systems 
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We calculate the resistivity of 2D electron (hole) gas, taking into account the degeneracy and the 
thermal correction due to the combined Peltier and Seebeck effects. The resistivity is found to be 
universal function of temperature, expressed in units of -^-(farZ) - . The giant parallel magnetire- 
sistivity found to result from the spin and, if exists, valley splitting of the energy spectrum. Our 
analysis of compressibility and thermopower points to thermodynamic nature of metal-insulator 
transition in 2D systems. 
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Recently, a great deal of interest has been focussed on 
the anomalous transport behavior of a wide variety of 
low-density 2D electron^, § and hole§, |, § systems. 
It has been found that, below some critical density, the 
cooling causes an increase in resistivity, whereas in the 
opposite high density case the resistivity decreases. An- 
other unusual property of dilute 2D systems is their enor- 
mous response to parallel magnetic field. At low temper- 
atures the magnetic field found to suppress the metallic 
behavior of 2D electron(hole) gas and result in strong in- 
creasing of resistivity upon enhancement of spin polariza- 
tion degree H f7|. At high temperatures the parallel mag- 
netoresistivity starts to be unaffected by magnetic field 
when the temperature exceeds a value being of the order 
of Zeeman energy. A strong perpendicular magnetic field, 
if applied simultaneously with the parallel one, results in 
suppression of parallel magnetoresistivityQ. Although 
numerous theories have been put forward to account for 
these effects, the origin of the above behavior is still the 
subject of a heated debate. 

The ohmic measurements known to be carried out at 
low current I — > in order to prevent heating. Usually, 
only the Joule heat is considered to be important. In con- 
trast to the Joule heat, the Peltier and Thomson effects 
are linear in current. As shown in Refs.|], [l(], 11 , the 
Peltier effect influences ohmic measurements and results 
in a correction to a measured resistance. When current 
is flowing, one of the sample contacts is heated, and the 
other cooled, because of the Peltier effect. The contact 
temperatures are different. The voltage drop across the 
circuit includes the Seebeck thermoelectromotive force, 
which is linear in current. Finally, there exists a thermal 
correction Ap, to the ohmic resistivity, p, of the sample. 
For degenerate electrons, Ap/p ~ (kT/p) 2 , where p is 
the Fermi energy. Hence, the correction may be compa- 
rable with the ohmic resistance of a sample when kT ~ p. 

In the present paper, we report on a study of low- 
T transport in 2D electron(hole) gas, taking into ac- 
count both the electron degeneracy and the Peltier-effect- 
induced correction to resistivity. |j [hJ . The parallel mag- 
netoresistivity found to originate from the spin and, if ex- 
ists, valley splitting of 2D energy spectrum. 

Let us consider, for clarity, the (100) MOSFET 2DEG 
system. The electrons are assumed to occupy the first 



quantum-well subband with isotropic energy spectrum 
e(k) = 4j^-, where m is the effective mass. The sample 
is connected (see Fig.^jA, inset) by means of two identical 
leads to the current source. Both contacts are ohmic. 
The voltage is measured between the open ends( "e" and 
"d" ) kept at the temperature of the external thermal 
reservoir. The sample is placed in a sample chamber 
with mean temperature To. 

According to our basic assumption, the contacts( "a" 
and " b" ) may have different respective temperatures T a 
and Tf,. Including the temperature gradient term, the 
current density j and the energy flux density q yield 

j = a(E-aVT), q = (aT - C/e)j - kVT, (1) 

where E = V£/e is the electric field, and ( — p — eip 
the electrochemical potential. Then, a = Ne 2 T/m is 
the conductivity, N the 2DEG density, k the thermal 
conductivity, and a the thermopower. 

In general, one can unambiguously solve Eq.(|I]), and, 
then find a difference of contact temperatures AT = 
T a —Tf, for an arbitrary circuit cooling. Since the electron- 
phonon coupling is weak below ~ IK, the heat con- 
duction from 2DEG to mixing chamber could predom- 
inately occur through the contacts of the sample and 
the leads connected to them. However, the experimen- 
tal observations demonstrate that 2DEG alone is the 
dominant thermal resistance in this problem. Actually, 
the cooling is provided by electron thermal conductivity 
found to follow Wiedemann-Franz law k = LTa, where 
L = is the Lorentz number. Accordingly, we further 
consider adiabatic cooling, with the 2DEG thermally in- 
sulated from the environment. Then, we will omit the 
Joule heating for actual I — ► case. We emphasize that 
under the above conditions, the sample is not heated. 
Indeed, at small currents, T a rs T, fz To. Hence, the 
amount of the Peltier heat, Q a — iAaTo, evolved at 
contact "a" and that absorbed at contact "b" are equal. 
Here, Aa is the difference of 2DEG and metal conduc- 
tor thermopowers. If it is recalled that the energy flux 
is continuous at each contact, the difference of the con- 
tact temperatures yields §|,JlO) AT = xf^T where Z n 



and w are respectively the sample length and width. Us- 
ing Eq.dfi), the voltage drop between ends "e" and "d" 
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is given by U = RI + AaAT, where R is the ohmic re- 
sistance of the circuit. The second term denotes conven- 
tional Seebeck thermoelectromotive force. Since AT ~ /, 
we finally obtain the total 2DEG resistivity as follows 



p(l+a 2 /L), 



(2) 



where p = 1/a is the 2DEG ohmic resistivity. In Eq. (g) 
we ignore conductors resistances and take into account 
that for the actual case of the metal leads Aa ~ —a. 
Eq.(|]) can also then be applied for 2D hole gas. 

Using Gibbs statistics and the above energy spectrum, 
the 2DEG density N = -(P) yields 



N = iVb£Fo(l/0, 



(3) 



where ft = -fcT£ln(l + ex P(^)) is the 2DEG ther- 
modynamic potential, and £ = kT/fi the dimension- 
less temperature. Then, N = Dp is the density of 
strongly degenerate 2DEG, D = ^ the 2DEG density of 
states(factor 2 accounts valley degeneracy) and F n {z) the 
Fermi integral. In Fig.|l|A, we plot the temperature de- 
pendence of the dimensionless concentration n = N/Nq 
given by Eq.(||). In the classical Maxwell-Boltzman limit 
(£ < 0, |£| <C 1), the 2D electron density is thermally 
activated, and, therefore, n = |£| exp(— 1/ |£|). In the 
case of strongly degenerate electrons (£ <C 1), we ob- 
tain n = 1 + £exp(— 1/£). Then, at elevated temper- 
atures £ > 1, the dependence of the 2DEG density 
n = 1/2 + £ In 2 becomes linear in temperature. It is 
to be noted that at finite temperature the 2DEG density 
always exceeds( see Fig.^, inset b) the zero-temperature 
value, i.e N > N . As a consequence, the density, de- 
duced from the period of Shubnikov-de-Haas quantum 
oscillations, Ng^^ = N , is less than that, Njjaii = N, 
obtained from low- field Hall resistivity data( see Ref . |l3|] ) 
in consistent with experiments. E3] We emphasize that 
the 2D electron density is a monotonic function of tem- 
perature (see Fig.|]A). Therefore, one might expect that 
the ohmic resistivity p(T) ~ 1 /N decreases with increas- 
ing temperature at constant carrier mobility. We now 
demonstrate that the total resistivity specified by Eq. (j^) 
can, nevertheless, increase in a certain temperature range 
owing to thermal correction. 

Following the conventional Boltzman equation formal- 
ism, the explicit formulae for the 2DEG thermopower 
(for the 3D case, see Pisarenko, 1940) can be written as 



2Ji(l/Q 1 
F (1/0 £ 



(4) 



Here, we assume, for simplicity, that the electron scat- 
tering is characterized by energy-independent momen- 
tum relaxation time. In the classical Maxwell-Boltzman 
limit (£ < 0, |£| <C 1) the thermopower is given by 
the conventional formulae a = — -(2 — l/£). Then, 
for strongly degenerated 2DEG (£ <§C 1), we obtain the 
temperature dependence of the thermopower (Fig.[l]B) 
as a = -|[7r 2 £/3 - (1 + 3£)exp(-l/£)]. At elevated 




FIG. 1: Zero-field 2DEG density(A) and thermopower(B) 
given by Eq.(^|) and Eq.(Q) respectively vs dimensionless tem- 
perature £. Asymptotes: |£ | -C 1 - dotted line, £ 3> 1 - dashed 
line. Low-temperature dependence of 2DEG density(a) and 
thermopower(b) given by Eq.(|^) at fixed spin polarization fac- 
tor B* = 0; 0.05; 0.25; 1.2. The arrows depict the magnetic 
field enhancement (diminution) of the thermopower(density). 
Inset(top panel): experimental setup. 

temperatures (£ > 1) the thermopower first grows with 
temperature, and then approaches the universal value 



k 2Fi(0) 
fo(0) 



- -| g^-j . Our support of the above be- 
havior is confirmed by low-T thermopower measurements 
datajD|, found to diverge at certain value ~ 0.6fc/e being 
of the order of a g ( see bold line in Fig.||, inset a). 

It's worth noting that the thermopower can be of 
the order of k/e. Accordingly, the thermal correction 
to resistivity may comparable with the ohmic resistiv- 
ity of 2DEG. In Fig.|2|, we plot the temperature depen- 
dence of the 2DEG resistivity given by Eq.(j^) at differ- 
ent Fermi temperatures Tp = p/k. At fixed tempera- 
ture, the resistivity increases with decreasing of 2DEG 
degeneracy. At certain Fermi energy (Tp — 0.25K in 
Fig.^) the T-dependence of the resistivity exhibits metal- 
lic behavior at T < Tp, and then becomes insulat- 
ing ( i.e. % < ) at T > T F . Within the low- 
temperature metallic region £ <C 1, the 2DEG resis- 
tivity can be approximated (see dashed line in Fig.H ) 
with p tot = p (l + tt 2 £ 2 /3), where p = ^(kpl)-^is 
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teracting 2DEG system Eq.@ yields #"(£) = DF^l/t), 
where F' (z) — rfF "( z ) is the derivative of the Fermi func- 
tion. At fixed temperature Fig]3| represents the depen- 



FIG. 2: Zero-field T-dependence of resistivity and ther- 
mopower(inset (a)), given by Eqs.(||^) for T F [A"]=2-0.25(step 
0.25), 0.2-0.05( step 0.05), 0.01, 0(bold line),-0.1,-0.2. Asymp- 
totes: { < 1 ■ dashed line, £ > 1 - dotted line for fixed 
Tf — 0.25K. Inset(b): 2DEG density vs Fermi energy at fixed 
temperature T — (piecewise bold line); 0.2; 0.4K. 



the resistivity at T — > 0, &f = y/2mfj,/h the Fermi 
vector, and I — kkpr/m the mean free path. Then, 
for the high-temperature (£ > 1) insulating region we 

obtain the asymptote p tot = Po 2+1/2 ' depicted in 
Fig.^J by dotted line. This results are confirmed by re- 
cent experiments flrj|, |f^| showing that for temperatures 
well below the Fermi temperature the metallic region 
data obey a scaling, where the disordered parameter kpl 
and dimcnsionless temperature T '/Tf appear explicitly. 
These experimental observations u& therefore rule out 
interactions, the shape of the potential well, spin-orbit 
effects and quantum interference effects |f7| as possible 
origins of the metallic behavior mechanism. It is to be 
noted that the conventional theory|l9| used to explain 2D 
metallic behavior|2(J is, however, failed to account both 
T — > and T > Tp cases. It is of particular interest the 
2DEG compressibility, K = = -jM, known to be a 
fundamental quantity generally more amenable to theo- 
retical and experimental analysis. ]2f|, 22, E3| For nonin- 
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FIG. 3: Dimensionless inverse compressibility vs Fermi tem- 
perature at zero magnetic field(thin lines) and hui c — lK(bold 
lines) at fixed temperatures T = 0.15, 0.25K. Dashed lines de- 
pict asymptotes at £ < 0, |£| <C 1 

dence of the inverse compressibility parameter, d(n) = 
e/Ke 2 , vs Fermi temperature. For strongly degenerated 
electrons (£ < 1) one obtains a constant value do = e/De 2 
in consistent with experimental findings. [B Then, upon 
decreasing of 2DEG density( i.e. [i — > 0) the experi- 
mental data^l], exhibit diminish and, further- 
more, the negative inverse compressibility compared to 
do. Usually, this behavior is explained ]22[ in terms of 
conventional Hartree-Fock exchange omitted in our sim- 
ple approach. However, for extremely low 2D density the 
inverse compressibility data always exhibits an abrupt 
upturn which cannot be explained^) within Hartree- 
Fock scenario. We argue that the above feature has the 
natural explanation within our model(see dashed line in 
Fig.|) since d = d exp{-l/ at £ < 0, |£| < I and, 
hence, exhibits T-activated behavior. 

In general case of 2DEG placed in strong perpendicular 
magnetic field B± the compressibility yields 
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K = 



D 



4 ^ „ cosh i 



2kT 



D 



(-l) k kcos(2irvk) 
sinh(27r 2 fc^) 



(5) 



Here, fi = -kTT £ In (l + cxp { t T ^ L )) is the ther- 

n 

modynamic potential modified with respect to 2DEG 
spin-unresolved Landau level(LL) energy spectrum s n — 
hu c (n + 1/2), where n — 0, 1.. is the LL number, oj c — 
the cyclotron frequency, and T the zero-width LL 
density of states. Then, v = -J*- is the filing factor. 
According to Eq.(||), at fixed magnetic field the depen- 
dence d(n) can be viewed( see Fig]3j) as a superposition of 
zero-field dependence and LLs related oscillations. Then, 
at £ < 0, |£| < 1 in presence of strong magnetic field 
Tilo c >> kT{ i.c.^v « 1) the inverse compressibility ex- 
hibits T-activated behavior as d = do^uexp — 

similar to zero- field case( see Fig.^j) . These findings are 
consistent with experimental observations. p2[ p3| 

We now address the question of parallel field mag- 
netoresistivity. The energy spectrum modified with re- 
spect to valley (24] A V [K] = A° + 0.6B||[T] and spin 
A a [if] = g*HBB\\ = 2.6S||[T] splitting is represented 
in Fig.fpa,, inset. For simplicity, we neglect further the 
density dependence of <?*-factor. Then, addressing the 
unresolved problem to extract zero-field valley splitting 
within low-field Shubnikov de-Haas measurements, we as- 
sume that A° = 0. Actually, these assumptions do not 
substantially affect our basic results. Note that the low- 
est spin-up states are always valley-split since the we have 
a ratio = A v /A s = 0.23 < 1. 

With the help of Gibbs approach, the explicit formulae 
for 2DEG density and thermopower yields 



N 



Not 



l-BU 



(6) 



2£^i(^)+ffb(^) 



l-Ei 



where E l = 0,B*,f3B*,(l + f3)B* is the energy deficit 
between the bottom of spin and(or) valley subbands and 
that of the ground state( see Fig^Ja, inset), B* = A s /fi 
the spin polarization factor. The system becomes com- 
pletely spin-polarized when B* = l(at T = 0). Let us dis- 
cuss qualitatively the 2DEG magnetoresitivity. At fixed 
temperature £ <C 1 the total number of electrons de- 
creases ( see Fig|j]a) upon increasing of the parallel mag- 
netic field as n = 1 - (1 + (3)B* /2 when B* < 1. Simul- 
taneously, the thermopower exhibits a strong enhance- 
ment see Fig.[j]b) when £ < B*. We argue, that both the 
enhancement of thermopower and diminish of 2DEG den- 
sity is responsible for the observed giant magnetoresistiv- 
ity in 2DEG(Fig.0b). In the opposite spin-polarized case 




A , K 



FIG. 4: (a)T-dependence of magnetoresistivity at fixed 
Tf = 0.25A" for different degree of spin polarization( from 
top to bottom) B* = 0.2; 0.15; 0.1 - 0.02(step 0.02); 0. 
Inset: energy diagram for valley and(or) spin subbands. 
(b)Magnetoresistivity vs Zeeman splitting at fixed temper- 
ature T — 0.07K for Fermi energies denoted in Fig.^. Dotted 
line represents the condition B* = 1 for spin-polarized 2DEG. 

B* > 1, the 2D electrons occupy the lowest valley-split 
subbands (Fig^,, inset), therefore n = 1/2 — /3£?*/4 ~ 
1/2. Then, for spin-polarized system the thermopower 
saturates( see Fig.|l|b), and, hence results in subsequent 
saturation of the parallel magnetoresistivity(see Fig.|]b). 
It worthwhile to mention that the strong parallel mag- 
netoresistivity observed in p-GaAs/AlGaAs system |2q ] 
can be explained within our model as well. However, the 
high-field behavior of resistivity remains unclear. 

At fixed magnetic field the T-dependence of magne- 
toresistivity given by Eqs.(||),(||) is similar to that rep- 
resented in Fig.|| for B=0 case. Then, for fixed Fermi 
temperature( see Fig.[|a) the above dependence becomes 
unaffected by applied magnetic field above certain tem- 
perature being of the order of Zeeman energy. This find- 
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ing is consistent with experimental observations. H It is 
worth to be noted that the conventional theory|ll| agrees 
with magnetoresistivity dataQ only qualitatively. 

Let us analyze in more detail the cooling conditions, 
which are known to influence the thermal correction to 
resistivity. ||, [To) It will be recalled that in the case of 
adiabatic cooling the electron temperature differs from 
the bath temperature To. We now consider the oppo- 
site situation of electron cooling due to, for example, 
finite strength of electron-phonon coupling. Following 
Ref.|(|, below -0.6K in Si MOSFET's the electron-to- 
phonon thermal exchange is given by a(T 3 — Tq ), where 
a =2.2 xlCT 8 W/K 3 cm 2 . When T - T < T , the ther- 
mal correction to the resistivity Ap is suppressed [|l0] by 
the factor 7 = ISliiiA^ where A = / ?ov / 3a/4K is a di- 
mensionless parameter. Actually, A is the ratio of out- 
going and internal heat fluxes associated with phonon 
related thermal leakage and electron heat diffusion, re- 
spectively. When A 1, the local cooling due to phonons 
can be neglected and the adiabatic approach is well jus- 
tified. In the opposite case of intensive cooling (A ^> 1), 
the difference of the contact temperatures AT becomes 
smaller, and, therefore, the thermal correction to re- 
sistivity Ap vanishes. For Iq — 1mm, To = 50mK, 
a = 2e 2 /h = 8 x 10 _5 Ohm _1 ( typical critical region 
conductance) we obtain A = 1.2, and, therefore, 7 = 0.8. 
It worthwhile to notice that Peltier correction to resistiv- 
ity becomes greater at ultra-low temperatures for short 
samples, since A ~ IqTq. 

We emphasize that both dc and ac ohmic measure- 



ments lead to a thermal correction. The correction is, 
however, strongly damped at high frequencies because of 
the thermal inertial effects. As demonstrated in Ref.||, 
the above quasi-static approach is valid below some crit- 
ical frequency f cr = x/Iq- For example, for degenerate 
electrons the thermal diffusion coefficient \ is of the order 

of the diffusion coefficient D = ^ ■ Assuming 

a = e 2 /h, Iq = 1mm, for GaAs-based structure we ob- 
tain x ~ Ti/m, hence f cr — 1.5kHz. We suggest that the 
spectral dependence of the 2D resistivity can be used to 
estimate the thermal correction. 

In conclusion, low-temperature ohmic measurements 
of a 2D electron (hole) gas involve a thermal correction 
caused by the Peltier effect. The resistivity of 2DEG 
with thermal correction included is found to be universal 
function of temperature, expressed in units h/e 2 (kpl)~ 1 ■ 
This universal behavior correlates with that found in ex- 
periments. Strong increasing, and, then subsequent satu- 
ration in parallel magnetiresistivity in Si-MOSFET found 
to results from both the spin and valley splitting of 2DEG 
energy spectrum compared with B = case. We suggest 
ac measurements as a powerful tool illuminating the im- 
portance of thermal correction to resistivity. Our analysis 
of compressibility and thcrmopower points to thermody- 
namic nature of metal-insulator transition in 2D systems. 
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